Boundary layers in rotating weakly turbulent Rayleigh-Benard convection. 
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The effect of rotation on the boundary layers (BLs) in a Rayleigh-Benard (RB) system at a 
relatively low Rayleigh number, i.e. Ra = 4 x 10 7 , is studied for different Pr by direct numerical 
simulations and the results are compared with laminar BL theory. In this regime we find a smooth 
onset of the heat transfer enhancement as function of increasing rotation rate. We study this regime 
in detail and introduce a model based on the Grossmann-Lohse theory to describe the heat transfer 
enhancement as function of the rotation rate for this relatively low Ra number regime and weak 
background rotation Ro > 1. The smooth onset of heat transfer enhancement observed here is in 
contrast to the sharp onset observed at larger Ra > 10 8 by Stevens et al. [Phys. Rev. Lett. 103, 
024503, 2009], although only a small shift in the Ra — Ro — Pr phase space is involved. 
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I. INTRODUCTION 



Normally the transition between different turbulent states is smooth, because the large random fluctuations that 
characterize the turbulent flow make sure that the entire phase space is explored and therefore the transitions between 
different states, that arc explored as a control parameter is changed, are washed out. A classical system to study 
turbulence is Rayleigh Benard (RB) convection [lj-Q- F° r given aspect ratio r = D/L (D is the cell diameter and 
L its height) and given geometry, its dynamics are determined by the Rayleigh number Ra = (3gAL 3 /(kv) and the 
Prandtl number Pr = v/k. Here /? is the thermal expansion coefficient, g the gravitational acceleration, A the 
temperature difference between the plates, and v and k are the kinematic and thermal diffusivity, respectively. The 
heat transfer in a RB system can satisfactory be described by the Grossmann-Lohse (GL) theory!]], and shows 
that RB convection has different turbulent regimes in the Ra — Pr phase space (see Fig. 3 of Ref. |l( ) . The case where 
the RB system is rotated around a vertical axis, i.e. Rotating Rayleigh-Benard (RRB) convection, at an angular 
speed fl is interesting for industrial applications and problems in geology, oceanography, climatology, and astronomy. 
The rotation rate of the system is non-dimensionalizcd in the form of the Rossby number Ro = \J f3gA/L/(2fl). The 
dynamics of RRB convection are thus determined by three control parameters, i.e. Ra, Pr, and Ro, and this leads to 
a huge Ra — Pr — Ro phase space, see Fig. [TJ 

It is widely understood § that rotation suppresses convective flow, and with it convective heat transport, when the 
rate of rotation is sufficiently large. However, experimental an d numerical [l6|-[24j studies on RRB convection 

have shown that rotation can also enhance the heat tr ansp ort with respect to the non-rotating case. This heat 
transport enhancement is caused b y E kman p umping [ill Il5l4l8l . I20L I22L I24T4261 ] and its efficiency depends strongly on 
the combination of Ra, Pr, and Ro |l7lll8l[24 |. In this paper we will discuss the results of Direct Numerical Simulations 
(DNS) that show that this heat transfer enhancement as function of the Ro number is smooth for relatively low Ra 
number, here Ra = 4 x 10 7 (see Fig. [^), while experimental and numerical data for Ra = 2.73 x 10 s and Pr = 6.26 
show a sharp onset for the heat transport enhancement, see Fig. ^jp. This difference is remarkable since only a small 
shift in the Ra — Pr — Ro phase space is involved (see Fig. [1]). 

In this paper we will first describe the flow characteristics found in the simulations. We will show that there is 
a smooth transition from one turbulent regime to another for the relatively low Ra number regime whereas a sharp 
transition is found for higher Ra. In section [TT] we will discuss the properties of the BLs found in the DNS in detail. 
Subsequently the laminar BL theory for flow over an infinitely large rotating disk will be discussed in section IIIII in 
order to explain the BL properties found in the DNS. The derived scaling laws from this theory will be used in a 
model based on the GL theory to describe the heat transfer enhancement as function of Ro for the relatively low Ra 
number regime with weak background rotation, see section ITVl 



II. NUMERICAL RESULTS FOR BOUNDARY LAYERS IN RRB CONVECTION 

The flow characteristics of RRB convection for Ra = 4x 10 7 , 1 < Ro < 00, and 0.2 < Pr < 20, are obtained from 
solving the three-dimensional Navicr-Stokcs equations within the Boussincsq approximation: 

Du ( Pr\ 1/2 , _ 1 _ 

7*=- VP+ U Vu + fe-^xu, (1) 



DO 1 



Dt (PrRa) 1 / 2 



V 2 6>, (2) 



with V • u = 0. Here z is the unit vector pointing in the opposite direction to gravity, D/Dt = dt + u • V the 
material derivative, u the velocity vector (with no-slip boundary conditions at all walls), and 9 the non-dimensional 
temperature, < 9 < 1. Finally, P is the reduced pressure (separated from its hydrostatic contribution, but containing 
the centripetal contributions): P = p — r 2 / (8R0 2 ), with r the distance to the rotation axis. The equations have been 
made non-dimensional by using, next to L and A, the free-fall velocity U = y//3gAL. A constant temperature 
boundary condition is used at the bottom and top plate and the side wall is adiabatic. Further details about the 
numerical procedure can be found in Rcfs. [27T - I291 ] . 

The first set of simulations is used to study the Ro number dependence of the following quantities: the normalized 
heat transfer, the thickness of the thermal BL, and the normalized averaged root mean square (rms) vertical velocity 
fluctuations. Here we have simulated RRB convection at several Ro numbers for three different Pr numbers (Pr = 0.7, 
Pr = 6.4, and Pr = 20). All these simulations are performed on a grid with 65 x 193 x 129 nodes, respectively, in 
the radial, azimuthal, and vertical directions, allowing for sufficient resolution in the bulk and the BL according the 
resolution criteria set in Ref. (30j . The Nusselt number is calculated in several ways as is discussed in detail in 
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Stevens etal. (2010 & 2009) 


• 


Liu and Ecke (2009 & 1997) 


Niemela et al. (2010) 


♦ 


Kunnenetal. (2008) 


Zhong and Ahlers (2009) 




Oresta et al. (2007) 


King et al. exp (2009) 


♦ 


Kunnenetal. (2006) 


King et al. sim (2009) 




Julien et al. (1996) 


Schmitz and Tilgner (2009) 


• 


Rossby (1969) 



FIG. 1. Phase diagram in Ra-Pr-Ro space for RRB convection. The data points indicate where Nu has been measured or 
numerically calculated. The data are obtained in a cylindrical cell with aspect ratio r = 1 with no slip boundary conditions, 
unless mentioned otherwise. The data from direct numerical simulations and experiments are indicated by diamonds and dots, 
respectively. The data sets, which are ordered chronologically in panel e, are from: Stevens et al. (2010 & 2009) fPH , Il8l . [24j] 
and the simulations of this study; Niemela et al. (2010) [H (F = 0.5); Zhong and Ahlers (2009) [jj.GJ]; King et al. (2009) 
"ill ; Schmitz and Tilgner (2009) [13] (free slip boundary conditions and horizontally periodic); Liu and Ecke (2009 & 1997) 
IliOl (square with T = 0.78); Kunnen et aL (2008) Oresta et al. (2007) (r = 0.5); Kunnen et al. (2006) [H (r = 2, 
horizontally periodic); Julien et al. (1996) [20] (F = 2, horizontally periodic); Rossby (1969) (varying aspect ratio). Panel a 
shows a three dimensional view on the phase space (see also the movie in the supplementary material), b) projection on the 
Ra-Pr phase space, c) projection on the Ra-Ro phase space, d) projection on the Pr-Ro phase space. 
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FIG. 2. Normalized heat transfer as a function of the Rossby number, a) Numerical results for Ra = 4 x 10 7 and Pr — 6.4. 
The open black squares and stars indicate the results obtained on the 193 x 65 x 129 and 257 x 97 x 193 grid, respectively, b) 
Ra — 2.73 x 10 s and Pr — 6.26. Red solid circles: experimental data Open black squares: numerical results |17j. Fl3|. The 
vertical dashed line indicates the position of the onset. 



Ref. [30| and its statistical convergence has been controlled. Some data for Pr = 6.4 have already been published in 
Ref. [18| . There the average was over 4000 dimensionlcss time units. The new results for Pr = 0.7 and Pr = 20 are 
averaged over 2500 dimensionlcss time units. Note that we simulated the flow for a large number of eddy turnover 
times to reduce the statistical error in the obtained Nusselt number results and to prevent the influence of transient 
effects. This is necessary to accurately resolve the transition regime where the heat transfer starts to increase and to 
accurately determine the flow statistics. Furthermore, we note that all simulations are started from a new flow field 
in order to rule out hysteresis effects. 

The second set of simulations is used to study the Pr number dependence of the same set of quantities. Here 
we simulated RRB for several Pr numbers and three different Ro numbers (Ro = oo. Ro ~ 3, and Ro = 1). The 
simulations for Pr ^ 0.70 are performed on a 97 x 257 x 193 and the simulations at Pr = 0.25 and Pr = 0.45 are 
performed on a 129 x 385 x 257 grid. The finer resolution is needed here as the structure of the flow changes and the 
Reynolds number based on the LSC increases for lower Pr. For most cases the flow is simulated for 400 dimensionless 
time units and 200 dimensionless time units were simulated before data are collected to prevent any influence of 
transient effects (30[. For Pr = 2 and Pr = 4.4 we averaged over 1200 dimensionless time units in order to obtain 
more accurate statistics on the velocity field. Note that the second set of simulations is partially overlapping with 
the first set of simulations. Wc find that results obtained on the different grids, are very similar, i.e. the difference is 
generally between 0.5% and 1%, see Figs. [2^,, [SJi and [7^,. Again all simulations are started from a new flow field in 
order to rule out any hysteresis effects. 

Fig. [3] shows the azimuthally averaged temperature profile at the cylinder axis for different system parameters. In 
previous numerical studies concerning (rotating) RB convection the thermal BL thickness is usually defined by either 
looking at the maximum rms value of the temperature fluctuations or by considering the BL thickness based on the 
slope of the mean temperature profile. In the latter case it is usually assumed that no mean temperature gradient 
exists in the bulk (the BL thickness according to this assumption is denoted by X S g~ ng ). The temperature gradient 
in the bulk is, however, strongly influenced by rotation [171 l33t| . and when rotation is present also by Pr, see Fig. 
2J We prefer to define the thermal BL thickness Xf as the intersection point between the linear extrapolation of the 
temperature gradient at the plate with the behaviour found in the bulk, see Fig. [3^. From now on this definition of 
the thermal BL thickness will be used here. 

For the relatively low Ra number regime, here Ra = 4 x 10 , the heat transfer enhancement as function of Ro is 
smooth, see Fig. [5^,. Note that although the behaviour of Xf (see fig. Eh), i.e. the horizontally averaged value of 
the radially dependent thermal BL thickness (Xf(r)), as function of Ro is similar for all Pr the behaviour of Nu is 
very different. This is due to the influence of Pr on the effect of Ekman pumping [Tt], [24j|. At low Pr the larger 
thermal diffusivity limits the effect of Ekman pumping and causes a larger destabilizing temperature gradient in the 
bulk [l7], Hil , see Fig. @] Due to the limited effect of Ekman pumping there is no heat transport enhancement for low 
Pr, see Fig. O 
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FIG. 3. a) Visualization of the definitions of the thermal BL thicknesses and X s e l ~ n9 (assuming zero gradient in the bulk), 
b-d) Azimuthally averaged temperature profiles for Ra — 4 x 10 7 and different Pr and Ro numbers at the cylinder axis r = 
for b) Ro = oo, c) Ro = 3, and d) Ro = 1. Black, blue, red, and dark green indicate the profiles for Pr = 0.25, Pr = 0.7, 
Pr = 6.4, and Pr = 20, respectively. The dots indicate the data points obtained from the simulations. The insets show the 
profile over the full cell. 



Fig. IHlshows that the volume averaged Re z ,rms (dimensionless rms velocity of the axial velocity fluctuations), which 
is a measure for the strength of the LSC, decreases strongly for strong enough rotation. The vertical dashed lines 
in Fig. [6] indicate the position where i?e ZjI . ms (f2)/i?e 2 rms (0) becomes smaller than 1, which we use to indicate the 
point at which the LSC strength starts to decrease. This value is determined by extrapolating the behaviour observed 
at low Ro numbers to reduce the effect of the uncertainty in single data points. In Ref. |18[ (see figure 3 of that 
paper) we also used this method to indicate the position of the onset of the heat transfer enhancement in the high 
Ra number regime. However, for this lower Ra number we do not find any evidence for a sudden onset around this 
point, see Fig. [5] where the vertical lines are plotted at the same positions as in Fig. [51 In contrast to the decrease in 
the volume averaged value of Re Ztrms the horizontally averaged value of Re z ^ rms at the edge of the thermal BL (thus 
at the distance Xf(r) from either the top or bottom boundary) increases. This indicates that Ekman pumping, which 
is responsible for the increase in Nu, sets in and no sign of Ekman pumping prior to the decrease in LSC strength is 
found. Fig. [5]thus shows that the flow makes a transition between two different turbulent states, i.e. a transition from 
a LSC dominated regime to an Ekman pumping dominated regime fl8j . Fig. [S^ shows no increase in the horizontally 
averaged value of Re z , rms at the edge of the thermal BL for Pr = 0.7, because the flow is suppressed for higher Ro, 
i.e. lower rotation rate, when Pr is lower, see the discussion in Ref. |24| . 
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FIG. 4. a) Horizontally averaged temperature gradient at midheight as a function of Ro for Ra = 4 x 10 7 and different Pr. Red 
open circles, blue diamonds, and black squares are the data for Pr = 0.7, Pr = 6.4, and Pr = 20, respectively, b) Temperature 
gradient at midheight as function of Pr for Ra = 4 x 10 7 and different Ro. Red open circles, blue diamonds, and black squares 
are the data for Ro — 1, Ro — 3, and Ro = oo, respectively. The vertical dashed lines in both graphs represent the point at 
which the LSC strength starts to decrease, see Fig. [6] 
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FIG. 5. a) Normalized heat transfer as a function of Ro for Ra — 4 x fO 7 and different Pr. Red open circles, blue diamonds, 
and black squares are the data for Pr — 0.7, Pr — 6.4, and Pr = 20, respectively. The dash-dotted line is the fit obtained 
by the model Eg. I221 with a — 55. The stars indicate the results from the grid refinement check (see text), b) Horizontally 
averaged thickness of the thermal BL \f as function of Ro. Symbols as in Fig. [5^, and the dash-dotted line indicates the 
same model fit as in Fig. [SJi. The vertical dashed lines in both graphs represent the point at which the LSC strength starts to 
decrease, see Fig. [6] 



The Pr number dependence of the Nu number and the thickness of the thermal BL is shown in Fig. [JJ From Fig. 
[7k, we can conclude that hardly any Pr number dependence on the Nu number exists in the weak rotating regime 
(Ro = 3). However, a strong Pr number effect appears for stronger rotation rates, where Ekman pumping is the 
dominant effect [i7j . Fig. [7}d shows that the effect of weak background rotation on the thermal BL thickness is largest 
for Pr f=s 2. The Pr number dependence of Re z ^ rms is shown in Fig. [HI The difference between the data points obtained 
for the bottom and top BL in Fig. [8] indicate the uncertainty in the results. Increasing the averaging time, which we 
checked for Pr = 2 and Pr = 4.4, reduces the differences for the data points obtained for the bottom and top BL. 
Furthermore, we note that it is important to take the radial thermal BL dependence (A|'(r)) into account for lower 
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FIG. 6. The normalized averaged rms vertical velocity fluctuations i?e Zjrms for Ra = 4 x 10 7 and different Pr as function of Ro. 
The black line indicates the volume averaged value of Re ZyTms . The red and the blue line indicate the horizontally averaged 
values of the Re z>rrns at a distance Xf(r) from the lower and upper plate, respectively. The vertical dashed lines again indicate 
the position where the LSC strength starts to decrease, a) Pr = 0.7, b) Pr = 6.4, c) Pr = 20. 
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FIG. 7. a) Normalized heat transfer as a function of Pr number for Ra = 4x 10 7 and different Ro. Red open circles and blue 
diamonds indicate the data for Ro — 1, and Ro = 3, respectively. The symbols (circles and diamonds) and stars indicate the 
results obtained on the 97 x 257 x 193 and the 65 x 193 x 129 grid, respectively, b) Horizontally averaged thermal BL thickness 
A| ! as function of Pr. Symbols as in a) and black squares for Ro = oo. 



Pr, where the radial BL dependence is strongest, and we excluded the region close to the sidewall (0.45 < r < 0.5) 
from the horizontal averaging in order to eliminate the effect of the sidewall. 

The computation of the kinetic BL thickness can be either based on the position of the maximum rms value of the 
azimuthal velocity fluctuations [lU, [34|, or on the position of the maximum value of e u := u • V 2 u, i.e. two times 
the height at which this quantity is highest, as shown in (30| . In ref. (35| we will compare in detail the profiles of e u 
and e u and will show that the latter is indeed suited to define the BL thickness. Here we first average e u horizontally 
in the range 0.05 < r < 0.45 before we determine the position of the maximum. This r range has been taken to 
exclude the region close to the sidewall, where e u misrepresents the kinetic BL thickness due to the rising plumes, 
and the region close to the cylinder axis, since there it is numerically very difficult to reliably calculate e' u , due to 

the singularity in the coordinate system. When the radially dependent kinetic BL thickness (A«"(r)) is horizontally 
averaged a small difference, depending on the averaging time, is observed between the bottom and top because of 
the specific orientation of the LSC. We note that the same quantity e u is used in Ref. J3(j, where it is shown that 
the kinetic BL thickness based on e u represents the BL thickness better than the one considering the maximum rms 
velocity fluctuations, which is normally used in the literature. The volume averaged value of e u is the same as the 
volume averaged kinetic energy dissipation rate e u (although it differs locally), which can easily be derived using 
Gauss's theorem [35[ . Fig. [9] shows the azimuthally averaged profiles for e u air~ 0.25L and in Fig. [10] the kinetic 
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FIG. 8. The normalized averaged rms vertical velocity fluctuations Re z , r ms for Ra = 4 x 10 7 and different Ro as function of 
Pr. The black line indicates the volume averaged value of Re z ,rms- The vertical error bars at Pr = 2 and Pr = 4.4 indicate 
the difference in the volume averaged value obtained after 1200 dimensionless time units (data point) and 400 dimensionless 
time units. The red and the blue line indicate the horizontally averaged Re z ,rms at a distance Xf(r) from the lower and upper 
plate, respectively, a) Ro = 3, b) Ro — 1. 
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FIG. 9. Azimuthally averaged profile of t„ for Ra = 4 x 10 7 and different Pr and Ro at the radial position r = 0.25L. The 
figures, from left to right, are for Ro = oo, Ro = 3, and Ro = 1. Black, blue, red, and green indicate the profiles for Pr = 0.25, 
Pr = 0.7, Pr — 6.4, and Pr = 20, respectively. The dots indicate the data points obtained from the simulations. The insets 
shows the profile over the whole domain. 



BL thickness based on the position of the maximum kinetic dissipation rate is shown as function of Ro and Pr. To 
compare the relative changes in the kinetic BL thicknesses the values are normalized by values for the non-rotating 
case. For all Pr numbers there is a change in the BL behaviour at the point where the LSC decreases in strength 
(vertical dashed lines in Fig. HUk). 

We conclude this section with a brief summary of the results obtained for the high Ra number regime [U LH| ■ For 
Ra > 1 x 10 8 a sudden onset at Ro = Ro c in the heat transport enhancement occurs, see Fig. [2b, where we find a 
smooth transition at lower Ra, see Fig. [5^. Fig. [TTk shows the volume averaged ratio Re Z!rms (Ct)/Re Zt rms(0)- The 
behaviour is similar to the one observed at lower Ra, see Fig. [51 For the high Ra number regime the onset at Ro c is 
defined as the point where the ratio Re Zjrms (fl) / Re z ^ rms (0) becomes smaller than 1. For the low Ra number this point 
indicates the position where the LSC strength starts to decrease and just as for the relatively low Ra number regime, 
Ekman pumping in the high Ra number regime is indicated by an increase of the horizontally averaged Re z _ rms value 
at the edge of the thermal BL. Although the two cases, i.e. the relatively low Ra number regime and the high Ra 
number regime both, show a transition between two different turbulent states the important difference between the 
two is that the transition is sharp in the high Ra number regime and smooth in the relatively low Ra number regime. 
The onset in the high Ra number regime is also observed in the behaviour of the BLs. Fig. [TTb shows that the kinetic 
BL thickness does not change below onset (Ro > Ro c ) and above onset the BL behaviour is dominated by rotational 
effects and thus Ekman scaling (proportional to Ro 1 ' 2 ) is observed. This scaling factor is also found in the laminar 
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FIG. 10. a) The thickness of the kinetic BL, average of bottom and top BL, based on the position of the maximum value of e u 
for Ra = 4 x 10 7 and different Pr. Red open circles, blue diamonds, and black squares are the data for Pr = 0.7, Pr — 6.4, and 
Pr = 20, respectively. The dash-dotted line indicates i?o 1//2 scaling, b) As in a), but now for different Ro. Red open circles, 
blue diamonds, and black squares are the data for Ro = 1, Ro = 3, and Ro = oo, respectively. The dash-dotted line indicates 
Pr 1//3 scaling. 



BL theory, which will be discussed in the next section. Finally, fig. [TTb shows the thermal BL thickness Xf. 
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FIG. 11. a) The averaged rms vertical velocities i?e z , rms as function of Ro. The black line indicates the volume averaged value 
of Re z ,rms. The red and the blue line, respectively, indicate the horizontally averaged Re z , r ms at the distance A|'(r) from the 
bottom and top plate, b) The thickness of the kinetic top and bottom BLs based on the position of the maximum e u , blue (red) 
for top (bottom) BL. The vertical dashed lines in panel a and b indicate the position of Ro c . Here the BL behaviour changes 
from Prandtl-Blasius (right) to Ekman (left) [l^|- c ) The blue circles (red diamonds) indicate \f for the bottom (top) plate. 



III. BOUNDARY LAYER THEORY FOR WEAK BACKGROUND ROTATION 

In the previous section we observed that there is a smooth increase in the heat transfer as function of the Ro number 
when the Ra number is relatively low. In this section we set out to account for the increase in the heat transfer as 
function of the rotation rate within a model, which extends the ideas of the GL-theory to the rotating case. In the GL 
theory the Prandtl-Blasius BL theory for laminar flow over an infinitely large plate was employed in order to estimate 
the thicknesses of the kinetic and thermal BLs, and the kinetic and thermal dissipation rates. These results were then 
connected with the Ra and Pr number dependence of the Nusselt number. In perfect analogy, in the present paper 
we apply laminar BL theory for the flow over an infinitely large rotating plate to study the effect of rotation on the 
scaling laws. We stress that employing the results of laminar BL theory over an infinite rotating disk to the rotating 
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RB case in a closed cylinder is fully analogous to employing Prandtl-Blasius BL theory for flow over an infinite plate 
to the standard RB case without rotation, where the method was very successful 0, 0-0] • 

In both cases the equations used to derive the scaling laws are time independent and therefore the resulting solutions 
are associated with laminar flow. However, evidently, high Rayleigh number thermal convection is time dependent. 
Therefore one wonders whether the derived scaling laws still hold for time-dependent flow over an infinite rotating 
disk. We will show that the Ro and Pr scaling that is derived is not changed when temporal changes are included. 
This is again in perfect analogy to the Prandtl-Blasius BL case where the scaling laws also hold for time dependent 
flow provided that the viscous BL does not break down 0. Indeed, recent experiments and numerical simulations 
[3a - t38| have shown that in non-rotating RB the BLs scaling wise behave as in laminar flow and therefore we feel 
confident to assume the same for the weakly rotating case. The basic idea of the model we introduce is to combine 
the effect of the LSC roll, which is implemented in the GL theory by the use of laminar Prandtl-Blasius BL theory 
over an infinitely large plate, and the influence of the rotation on the thermal BL. 

The system we are analyzing to study the influence of rotation on the thermal BL thickness above a heated plate 
is schematically shown in Fig. 112b . It is the laminar flow of fluid over an infinite rotating disk. The disk rotates with 
an angular velocity Qd and the fluid at infinity with angular velocity Qp = sflo, with s < 1. Fig. 112b shows that a 
positive radial velocity is created due to the action of the centrifugal force. Because of continuity there is a negative 
axial velocity, i.e. fluid is flowing towards the disk. The system has been analyzed before in the literature, e.g. Rcfs. 
[39l - l43j . Here we will briefly summarize the procedure. 

The system is analyzed by using the Navicr-Stokcs equations in cylindrical coordinates and assuming a steady 
stationary, axial symmetric solution. To reduce the Navier-Stokes equations to a set of ordinary differential equations 
(ODE) we employ self-similarity. The first step is to determine the dimensionless height in the system just as in the 
Prandtl-Blasius approach, in which the BL thickness scales as [44[ S ~ \fuxJ\J . For the case of a large rotating disk 
in a fluid rotating around an axis perpendicular to the d isk, the thickness of the BL can be estimated by replacing 
U by flox (44[. The thickness then scales as S ~ y/u/Qn- The similarity variable in the system is the dimensionless 
height 




According to von Karman, the following self-similarity ansatz for the velocity components and the pressure can be 
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taken mU 



u = rn D F(C), (4) 
v = rQ D G(0, (5) 



w = v^ff(C), (6) 
p = P M D P(0 + lps 2 n 2 D r 2 . (7) 

After substitution into the Navier-Stokes equations one obtains a system of four coupled ODEs, 

F 2 +FH-G 2 -F" + s 2 = 0, (8) 

2FG + HG - G" = 0, (9) 

P + HH - H" = 0, (10) 

2F + H =0, (11) 

where the prime indicates differentiation with respect to £. This set of ODEs must be supplemented by the boundary 
conditions 

u = 0, v = rfl D , w = forz = 0, (12) 
u = 0, v = srflD for z = oo. (13) 

When substituting the self similarity ansatz (01 - ([7]) into these boundary conditions one obtains 

F = 0, G = 1, H = for C = 0, (14) 
F = 0, G = s for C = oo. (15) 

Note that the boundary condition at infinity together with the continuity equation (|11[) gives H (£ — » oo) = 0. 
One can further simplify the set of ODEs by realizing that the ODE for the pressure (fT0|) is decoupled from the 
ODEs determining the velocity profiles by using the continuity equation (|11[) in (|10[) and subsequently integrating 
this relation. The velocity profiles, for the von Karman case (s = 0), are shown in Fig. IT2"b . The inset shows that 
the dimensionless kinetic BL thickness X u = X u /6 decreases with increasing relative rotation rate s of the fluid at 

—99% 

infinity. This is due to the decreasing effect of the centrifugal force. We determined X u , the dimensionless kinetic 

BL thickness at which the velocity has achieved 99% of the outer flow velocity, using the tangential velocity profile, 

s i 

i.e. when G(Q = s + 0.01(1 — s). Additionally, we calculated X u , the dimensionless kinetic BL thickness based on the 
slope of the tangential velocity at the disk. The scaling of the kinetic BL predicted by the above rotating BL theory, 
i.e. Ro 1 / 2 , is the classical Ekman scaling. In the simulations of RRB we find the same scaling of the kinetic BL once 
the flow is dominated by rotational effect, i.e. Ro < Ro c , see Fig. [TTb . 

The GL theory heavily builds on laminar Prandtl-Blasius BL theory, which describes the laminar flow over an 
infinite plate. In the Prandtl-Blasius theory the temperature field is assumed to be passive to derive the Pr number 
scaling. As we want to extend the GL theory to the rotating case we keep this analysis analogous to the Prandtl- 
Blasius theory. Therefore, we assume the temperature field to be passive in order to derive the scaling laws as function 
of the Pr number. We non-dimensionalize the temperature by 

0(C) = (16) 

where 6b is the temperature of the bottom disk, and Oao < 9b is the ambient temperature. Then one obtains the 
following ODE describing the temperature field j45l - l47| 

9" = PrH(Q9', (17) 

where the prime indicates a differentiation with respect to £. The boundary conditions are 

<9 = lforC = 0, (18) 
(? = 0forC = oo. (19) 
The resulting system of ODEs subjected to the boundary conditions, is solved numerically with a fourth order Runge- 
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Log(Pr) Log(Pr) 

(a) (b) 

FIG. 13. a) The thermal BL thickness Xg = X$/8 as function of Pr for the von Karman case (s = 0) in blue, for a — 0.5 in 
red, and for the Ekman case (s « l,Ro* = 10 -3 , see definition Eq. (|A2[l ) in black. The solid lines are for Xg and the dashed 

— 99% _ i /-a 

lines for Xg . Note that the scaling for the thermal BL thickness goes asymptotically to Pr ' in the high Pr regime and to 
Pr^ 1 in the low Pr regime. The inset shows Pr croBS , the transition between the high and the low Pr regime as function of s 
based on the behaviour Xg (solid line) and Xg 9% (dashed line). Note that the low Pr regime is more favored for higher s. b) 
The effective scaling exponent 7 in Xg ~ Pr 1 as function of Pr for the Blasius case (no rotation, green), the von Karman case 
(s = 0, blue), and the Ekman case (s « 1) for two Rossby numbers where black indicates the solution for Ro* = — 10~ 3 and 

s ; 99% 

magenta for Ro* = — 10 . The solid lines indicate the effective scaling for Xg , and the dotted lines for Xg . The red line 
indicates the analytic prediction, see Eq. (|A12[) . for the von Karman case. 



Kutta method using a Newton-Raphson root finding method to find the initial conditions. One can take the analytic 
solution for the Ekman case (s ~ 1), see appendix [XJ or the known solution for the Bodewadt case, see for example 
[111 Hi), as one of the starting cases to determine the solutions over the whole parameter range in s and Pr. 
In this way we obtain the full temperature profile for all s and Pr. For the heat transfer the most relevant quantity 

is the thermal BL thickness Xg, which we non-dimensionalizcd by S, thus Xe = Xe/S. One can distinguish between X e 

and X s ° , the dimcnsionless BL thickness based on the 99% criterion, thus when 6*(£) = 0.01. Fig. IT3a shows that the 

99% — s i 

asymptotic scaling of X s and Xg is the same. Furthermore, the figure shows that rotation does not influence the 
scaling of the thermal BL thickness in the high Pr regime, because the same scaling, namely proportional to Pr -1 / 3 , 
is found as for the Prandtl-Blasius case. However, the rotation docs influence the scaling in the low Pr regime, where 

now Xg ex. Pr^ 1 instead of Xg cx Pr^ 1 / 2 as found in the Prandtl-Blasius case. Notice that Xg ° > Ag , which is due to 

the decreasing temperature gradient with increasing height. 

In Fig. [T3"b we show the effective power-law exponent 7 = (dlog Xg)/ (d log Pr) of an assumed effective power law 

Xg ~ Pr 1 . It confirms that the effective scaling in the high Pr regime is the same for the Prandtl-Blasius (no rotation) 

and the von Karman case (s = 0), but already at Pr = la significant difference is observed. 

The temperature advection equation (( 1 T|) directly suggests the following relation between the scaling of the thermal 

s / s [ 

BL thickness X g and the scaling of the axial velocity at the edge of the thermal BL, Hbl ~ l/PrX g . This immediately 

implies for the low Pr regime, with Xg ~ Pr^ 1 , that Hbl is independent of Pr. For the high Pr regime, with 

Xg ~ Pr- 1 / 3 , it gives Hbl{Pt) ~ Pr 2 / 3 . The scaling of the thermal BL thickness in the low Pr regime can be 
understood on physical grounds. In this regime Xg 3> A u and, the kinetic BL is fully submerged in the thermal 
BL. The axial velocity at the edge of the kinetic BL (fT"5|) is Hbl = H{( — > 00), as can be shown by applying mass 
conservation expressed by Eq. (|11[) . As a consequence, in the low Pr regime the axial velocity is constant in almost 

the whole thermal BL. Then Eq. (|17[) can trivially be integrated and immediately gives X g ~ Pr^ 1 . This derivation 
is valid for all s, i.e. the scaling in the low Pr regime does not depend on the rotation of the fluid at infinity. The 
scaling in the high Pr regime is also independent of the rotation of the fluid at infinity. In the Ekman case s«l, see 
appendix lAl this 1/3 scaling regime shifts towards very large Pr. 

The equations (|5))- (fTTj) are time independent and therefore the resulting solutions are understood to describe laminar 
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flow. Temporal changes can easily be included by adding d^ug, dj:u r , and d^u z where t = tflo without changing the 
Ro and Pr scaling discussed above, i.e. the derived scaling laws above still hold for time dependent flow provided 
that the viscous BL does not break down. This is in perfect analogy to the Prandtl-Blasius BL case where the scaling 
laws also hold for time dependent flow Q. 

To investigate the crossover between the high and the low Pr regime we define Pr cross as the crossover point. 
Pr C ross is calculated by determining the intersection between the asymptotic behaviour of the high and the low Pr 

s [ 99% 

regime. To calculate Pr cross we considered X g (Pr) and X g (Pr). The inset of Fig. [T3"a shows that the low Pr regime 
becomes more favored when the rotation of the fluid at infinity becomes stronger. Then the kinetic BL thickness 
becomes thinner and the axial velocity decreases, i.e. the thermal BL thickness increases. Because the thermal BL 
thickness decreases with increasing Pr the crossover shifts towards higher Pr as is shown in the inset of Fig. [T3k . 
This effect is visible in RRB as shown in Fig. 5b of Ref. JUJ ■ Here it is shown that for the non- rotating case X u < X g l 
when Pr < 1. When the rotation rate is increased, i.e. Ro is lowered, this transition shifts towards higher Pr and for 
Ro = 0.1 A u < \g l when Pr < 9. 

In summary, the laminar rotating BL theory explains the Ro 1 / 2 scaling of the kinetic BL thickness in RRB convection 
and the shift of the position where \ u = Xg towards higher Pr when the flow is dominated by rotational effects. 



IV. MODEL FOR SMOOTH ONSET IN RRB CONVECTION 



In this section we will introduce a model in the spirit of the GL approach in order to describe the smooth increase 
in the heat transfer as a function of Ro that is observed for relatively low Ra and weak background rotation. Since the 
GL theory assumes smooth transitions between different turbulent states the model is limited to the relatively low Ra 
number regime, since a sharp onset as function of Ro is found for Ra > 1 x 10 8 . Furthermore, we assume that then 
the LSC, a basic ingredient of the GL model, is still present. In this simple model we neglect the influence of Ekman 
pumping, because it is a local effect that is rather insignificant at weak background rotation. This is supported by 
the results in Fig. [6] and the EPAPS document of Ref. [18j . where we find no evidence for Ekman pumping at weak 
background rotation. When strong rotation is applied Ekman pumping is the dominant effect and the validity regime 
of our model is thus restricted to weak background rotation. The basic idea of the model is to combine the effect 
of the LSC roll, which is implemented in the GL theory by the use of the laminar Prandtl-Blasius theory over an 
infinitely large plate, and the influence of rotation on the thermal BL. 

Applying laminar BL theory requires that the viscous BL above the flat rotating plate does not brake down. 
This assumption will now be verified. The stability for the von Karman flow has been studied theoretically and 
experimentally by Lingwood [48l - l50j . showing that instability occurs at Re w 510 where Re is defined as Re = r^J^ljv 
(with r the distance to the rotation axis). Lingwood also pointed out that other ex per imental studies show the same 
transition point within a very narrow Reynolds number range Re = 513 ± 15, see [48j and references therein. More 
recent experiments show similar results |5lll5^. For the Ekman case the instability occurs at Re ~ 200, see [Hoj . For 
the case under consideration it can be shown that Re < 55 (with r = 12.5 cm, £1 rj 0.20 rad/s, and v = 1 X 10 6 m 2 /s) 
[l7| . It can safely be conjectured that laminar BL theory can be applied as the estimated Reynolds number is an 
order of magnitude smaller than the critical Reynolds number. For further discussion on the stability of the rotational 
flow we refer to the classic Refs. [53j] . 

We introduce Xg R = X g (Pr)y' 'v /f2, see X g (Pr) in Fig. IT3k . as the thermal BL thickness based on the background 
rotation and Xgc as the BL thickness based on the LSC roll. Furthermore, T is the diameter-to-height aspect ratio of 
the RB cell and we set the radial length r = (TL)/2. Note that this is analogous to the length L which is introduced 
in the GL theory for the length of the plate. Thus the Reynolds number based on the background rotation is 



QL 2 T 2 1 , . 

Re « = — *R-o- (20) 

To calculate Xg R we used v = 1 x 10~ 6 to 2 /s (water) and T = 1 and we set Xgn = Xgc at Ro = oo. The strength of 
the LSC roll is taken constant and Xgc is known from Xgc/L = (2Nu)^ 1 . 

We model the increase of Nu as a crossover between a convection role dominated BL and a rotation dominated BL. 
Thus without rotation the BL thickness is determined by the LSC roll, i.e. Xg = Xgc, and when rotation becomes 
dominant the BL thickness is determined by the rotating BL thickness, i.e. Xg = Xg R . We now model the crossover 
between these two limiting cases as 

Xg_ _ \/Re R Xg R + a*\/Re c Xgc 
L ~ (^Re~^+a*^/Re^)L 
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Here the square root of Reynolds has been chosen since the dimensionlcss BL thicknesses scale with 1/y/Re. Wc 
rewrite the above equation in terms of Ro, using Ren oc Rec / Ro, 

\h — ts=Mr + Ae» c 

/x v ay Bo fOO\ 

and we use the free parameter a to fit the model with the numerical data shown in Fig. [5] It can be concluded 
that the presented model, based on the approach of the GL theory, indeed reflects the increase of Nu (as compared 
to the case without rotation) observed at relatively low Ra number (Ra = 4 x 10 7 ) when the LSC is still present. 
Furthermore, the thickness of the thermal BLs is also reflected correctly by the model. Note that the large value of 
the parameter a = 55 indicates that the influence of the rotation is rather weak before Ekman pumping sets in and it 
also explains that for higher Ra (Ra > 1 x 10 s ) no heat transfer enhancement is observed below onset. The sudden 
(instead of smooth) transition is then fully determined by the rotation rate where Ekman pumping sets in. This may 
be because at higher Ra the thermal BL is already much thinner due to the stronger LSC and therefore the effect 
of weak rotation is not sufficient to result in a significant thinner thermal BL. When Ro < Ro c the model cannot be 
used, since Ekman pumping is dominant in this regime which is responsible for the strong increase observed in Nu 
when Ro < Ro c . 



V. CONCLUSIONS 



To summarize, we have studied the effect of rotation on the RB system at relatively low Ra number, i.e. Ra = 4 x 10 7 
by using DNS. We find a smooth increase of the heat transfer as function of the rotation rate when weak rotation 
is applied. To describe this heat transfer enhancement we have extended the GL theory to the rotating case by 
studying the influence of rotation on the scaling of the thermal BL thickness. It is based on a similar approach as 
in the laminar Prantl-Blasius BL theory over an infinitely large plate, as we analyzed the flow over an infinitely 
large rotating disk where the fluid at infinity is allowed to rotate. Just as in the Prantdl-Blasius BL theory we used 
a passive temperature field to calculate the characteristics of the thermal BL. It turns out that weak background 
rotation does not influence the scaling of the BL thickness in the high Pr regime, because again Pr" 1 / 3 scaling is 
found. However, rotation does influence the scaling in the low Pr regime where we find a scaling of Pr" 1 instead of 
p r -i/2 f ounc j m the Prandtl-Blasius BL theory. With our model for the thermal BL thickness, see Eq. (|^|) . we can 
explain the increased heat transfer observed in the relatively low Ra number regime before the strength of the LSC 
decreases. The model neglects the effect of Ekman pumping as this effect is rather insignificant before the strength of 
the LSC decreases, i.e. the regime to which the model is applied. This means that the model cannot predict the heat 
transfer enhancement that is observed at moderate rotation rates where Ekman pumping is the dominant mechanism. 
The contrast between the smooth onset at Ra = 4 x 10 7 and the sharp onset at Ra > 1 x 10 8 is remarkable since 
only a small shift in the Ra — Pr — Ro phase space is involved. 

Acknowledgments: Wc thank R. Vcrzicco for providing us with the numerical code and F. Fontenele Araujo, F. 
Busse, G.J.F. van Hcijst, C. Sun, and L. van Wijngaardcn for discussions. The work is sponsored by the Foundation 
for Fundamental Research on Matter (FOM) and the National Computing Facilities (NCF), both sponsored by NWO. 
The numerical calculations have been performed on the Huygens cluster of SARA in Amsterdam. 



Appendix A: Ekman boundary layer theory 

In this appendix the results obtained from the model with weak background rotation will be compared with analytic 



results obtained from Ekman BL theory |54| . which uses a rotating reference frame. In the Ekman case the fluid at 
infinity is rotating at almost the same velocity as the disk, i.e. the limiting case of the model will be checked. Wc 
will indicate all quantities calculated in the rotatin g re ference frame with an asterisk. 

We will use the BEK model, presented in [50l [55ll56|. to derive a similar ODE as in section HTT1 for the temperature 
advection equation in the rotating reference frame. In the BEK model the following self-similarity assumption for the 
axial velocity is proposed: 



to = VM*Ro*H*((), 



(Al) 
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where £ = zy/Sl* /v. Here, n* is a system rotation rate, and Ro* is a constant determined by the rotation rate. We 
call this still the Rossby number since it also represents a dimcnsionlcss inverse rotation. In the BEK model Ro*, 
Af2, and 0*, respectively, are defined as 

Afl 

Ro* - W , (A2) 

An = n F - n D , (A3) 
w = -%- + -^-. (A4) 

2-Ro* 2 + Ro* y ' 

We obtain the following temperature advection equation in the BEK model 

9" = PrH*{i)Ro*9' . (A5) 
From Eqs. (|A2[) - (|A4[) one obtains, after some algebra, 

'2 + Ro* - Ro* 2 ' 



_ n F 
s ~ n^ 



2 - Ro* - Ro* 2 



(A6) 



This means for the Ekman case (s w 1, i.e. Ro* — > 0) that n F ~ nu(l + Ro*). Thus when the fluid at infinity is 
rotating slower than the disk Ro* is negative. From now on we assume Ro* to be negative, i.e. s < 1. 

Using Ekman BL theory (54j one can derive analytic solutions for the radial and tangential velocity profiles in the 



rotating reference frame. These analytical solutions read [541 ]: 



u E = -Afire^sinC = r n*Ro*F*(() , (A7) 
v E = Afir(l - e~ c cosC) = rn*Ro*G*(() , (A8) 

with C as in With the analytic expression for the radial velocity (|A7[) and the continuity equation one obtains 

w E = An. (l - e~ c [sinC + cosC]) 

V 

= VM*Ro*H*(C) . (A9) 

In the case Ro* — ^ it is found that n d m 0*, thus £ ~ C- In particular, the expression for the axial flow reduces to 
H*(£) = y/Q* /Qd(1 — [sin£ + cos£]), where yO* /Qd ~ 1- We find that the analytic expressions and the above 
numerical solutions are identical within numerical accuracy for the limiting case (Ekman solution, s ~ l,i?o* — > 0). 
(Note that a coordinate transformation has to be applied as the Ekman solution is expressed in the corotating reference 
frame, whereas the numerical solution has been defined in the laboratory frame.) 

Now we use this approach to calculate the BL characteristics for the Ekman layer. With the temperature advection 
equation (|A5|) and the analytic expression H*(£) for the axial velocity we determine the effective scaling exponent 7 
in Xg ~ Pr 1 as function of Pr for Ro* = — 10~ 3 and Ro* = — 10~ 5 . Fig. H3l shows that when Ro* goes to zero the 
low Pr regime (Xg 3> A u ) is extended to higher Pr, because the kinetic BL thickness becomes thinner and the axial 
velocity becomes smaller, i.e. the thermal BL becomes thicker. 

Substitution of H* (£) in (IA5|) yields the following expression for the temperature gradient in the Ekman layer 

— =Cl eti e ~ ic ° s ^) PrR °*\ (A10) 

where C\ is a constant of integration which does not depend on £. To determine the constant C\, and thereby Nu, one 
needs to integrate Eq. (|A10[) . An analytic result can be derived by substituting cos£ = (e 1 ^ + e~^)/2, A = PrRo*, 
z = B = (i — 1)/A (and B the complex conjugate of B), and evaluate the integral 

0(z) = 9l I e^^e^'dz + C 2 . (All) 



16 



The integration constants C\ and C2 arc determined by the boundary conditions 0(£ = 0) = 1 and 0(£ — > 00) = 0: 
C2 = and, for small A, 

G~*-k A -h A ' +0 <f)- (A12) 

The thermal BL thickness scales according to Xg oc e~ A /C\. With (|A12[) we immediately see that in the small 
Prandtl number limit: Xg rx Pr . Comparison with Fig. 113b reveals that the analytic results are in good agreement 
with the numerical data for Xg represented by the solid lines. Moreover, it predicts the scaling for Ro* = — 1 (the von 
Karman case, thus far outside the regime of applicability of the Ekman analysis) surprisingly well, see Fig. 113b . 
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